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Abstract 



' ■ In this letter new, closed and compact analytic expressions for the evaluation 

Q ■ of resonant energies, resonant bound-states, eigenvalues and eigenfunctions 



for both scattering and bounded n-cell systems are reported. It is shown 



^ . that for (scattering and bounded) 1-D systems the eigenfunctions ^^^(^(z) are 

^^ . simple and well defined functions of the Chebyshev polynomials of the second 

f~^ . kind Un, and the energy eigenvalues i?^ jy (in the /i-th band) are determined 

o. 

^D . by the zeros of these polynomials. New insights on the energy gap and the 



localization effect induced by phase coherence are shown. 
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Different techniques and approaches have been developed to determine, as exactly as 
possible, eigenvalues and eigenfunctions. In quantum theory the knowledge of these quan- 
tities represents a central aim and afford the possibility of evaluating other quantities. It is 
well known that, except for a small number of cases, it is not easy to obtain eigenvalues with 
the precision that one would like to have and to deduce explicit analytical expressions for the 
eigenfunctions. In the actual solid state physics theory, the eigenfunctions and eigenvalues 
are in some sense beyond it, and the theory is basically designed to describe systems in the 
continuous spectrum limit. The famous and widely accepted Bloch function e*'^^''-u^,fc(r), 
strictly valid only for infinite periodic systems[l], becomes the obvious and natural start- 
ing point for almost any approach of periodic systems, even though the periodic function 
M^fcg(r) remains in general practically unknown. Calculations based on the Bloch function 
lead also, almost automatically, to introduce the reciprocal space that inexorably, aside from 
text books' didactic presentations, may obscure the analysis. As a consequence, great efforts 
and imagination had to be bestown to explain the vast world of mechanisms, effects and 
properties discovered and populating the actual solid state and condensed matter physics[2- 
6]. The allowed and forbidden energy bands, which should not be confused with the energy 
eigenvalues of finite systems, are usually determined after a complex and perhaps painful 
numerical calculations with the aid of experimental input to fit parameters. On the other 
hand, high precision experiments and energy dependent applications, as those related to 
optical excitations, require one to solve the fine structure in the bands. This has been so far 
difficult for the actual theory. It is the purpose of this letter to extend the analytic methods 
of the scattering theory, for finite periodic systems with arbitrary potential shape, to deter- 
mine compact and general expressions for an easy and simple evaluation of the eigenvalues 
and eigenfunctions, which are neither the energy bands nor the Bloch functions. In other 
words, the present approach allows one to solve the fine structures in bands completely. 

Concerning band structures and electronic properties, there is an abundant literature 
and great contributions extending from the nearly free electron model to pseudopotential 
methods and Local Density Approximations. Important attempts were made forty years ago 



by Kohn to study the analytic properties of Bloch waves and Wannier functions[2]. Since 
those years, the complexity of involved Green's function and matching methods[7], together 
with a complete zoo of effective masses, perturbative methods, envelope functions, dispersion 
relations, and much more, have dominated the theory of real periodic systems [8- 11]. 

A completely different and natural approach to deal with finite (or locally) periodic 
systems is being developed[12]. In this approach the Bloch functions, reciprocal spaces, 
Brillouin zones, and other intricate concepts and knotty theoretical structures are not needed 
at all in order to calculate fundamental quantities. The theory of finite periodic systems 
(TFPS), based on the transfer matrix method (originally proposed to calculate scattering 
amplitudes, and related properties in superlattices), will be expanded to include also the 
stationary properties and to deduce, new and closed formulas for quantities as fundamental 
and basic as the eigenvalues and eigenf unctions. 

To encompass most of the various types of systems, I shall consider three representative 
cases, distinguished by their boundary conditions. Hence, I shall refer to open, bounded 
and quasi-bounded systems, as shown in figure 1. While for scattering systems (figure la)) 
I will refer to resonant energies and resonant bound-states, for bounded systems I will talk 
of energy eigenvalues and eigenfunctions. For the reasons of simplicity and lack of space, 
I will restrict to the widely used one channel (1-D) approximation. Some multichannel 
generalizations are almost straightforwardly obtained[13]. 

It is well known, in the scattering theory, that resonant transmission occurs precisely 
when the incident energy coincides with a bound-state energy in the scatterer system, which, 
for our purpose, is locally periodic with n cells (of length l^ each) in the transmission 
direction. For this kind of systems, the transmission amplitudes of iV-propagating modes 
(A^-channels) are obtained from[12] 

^^•" ^ PN,n-p-'aPpr,,n-i ^^^ 

where the functions pAr„ are N x N matrix polynomials fully determined in terms of the 
single-cell transfer matrix M, which for time reversal invariant systems has the structure 



M 



( a P^ 



(2) 



\I3* a* J 

In the one channel (one propagating mode) approximation a and (3 are complex scalars and 
Pi^n is the Chebyshev polynomial of the second kind Un{o.R) evaluated at the real part of a. 
Because of the close relation between the resonant structure and the spectral properties, we 
can determine simple expressions to evaluate the resonant energies. If we use the identity 
UnUn-2 = U^-i — 1) the whole n-cell system transmission coefficient can be rewritten as 

|f|2_ H /ON 

Here |t| is the single-cell transmission coefficient. It is clear from this expression that the 
transmission resonances occur when the polynomial f/„-i becomes zero. Therefore, the z/-th 
resonant energy E^j^^iy is a solution of 

(ttij)j. = cos — (4) 

f 1/ 

where (0:^)1^ is the z/-th zero of the Chebyshev polynomial and u = 1,2, ...n — 1. The 
index fi labels the bands, peculiar to periodic systems and entirely determined by the phase 
coherence. In the transfer matrix approach the allowed energy bands are those energies 
which satisfy the condition |a/?| < 1. To illustrate the use of equation (4), let us consider 
the square-barrier superlattice of figure la). This is similar to the familiar Kronig- Penney 
model (except that here we have a finite number of cells n) with barrier height Vo, and valley 
and barrier widths a and b, respectively. In this case 

COS ki,a cosh qi,b — sin ki,a sinh q^b = cos — (5) 

2Kq^ n 

with kl = 2nrlE^y/h^ and q^ = 2ml{Vo — E^^y)/h^ . Notice that each energy band contains 
the same number of resonant energies as confining wells have the scatterer periodic system. 
In this case is n — 1. In figure 2, some of these energies and the associated level densities 
p{E) are plotted for different values of n. The level density behavior tends rapidly, as a 



function of n, to that of the Kronig-Penney model, though the continuous spectrum hmit is 
reached only when n ^ oo. 

An important extension of the scattering approach, and the transfer matrix method, 
is to studying stationary properties of bounded periodic systems. If we have systems like 
those shown in figures Ib-c), we can apply the main results of the TFPS and introduce the 
boundary conditions. For a periodic n-cell system (with length nlc) bounded by infinite 
hard walls, it is easy to show that the energy eigenvalues are determined from 

(a„ -a: + P:-Pn) = (6) 

which, using the relations[12] «„ = f/„ — Q;*f/„_i and /9„ = /3f/„_ican be written as f/n-i(a/ — 
f3j) = 0. Here the subscript I refers to the imaginary part. Notice that n — 1 of the energy 
eigenvalues of the bounded system are the zeros of the Chebyshev polynomial f/„-i. For the 
particular example shown in figure lb), which length is nlc + a, the eigenvalues are obtained 
from UnSinka + (a/cos/ca — aRsinka — Pi)Un-i = 0. These eigenvalues are used below 
when evaluating the eigenfunctions. 

Another type of system that I would also like to consider is a periodic n-cell system 
confined by finite potential walls as shown in figure Ic). In this case, assuming E < V^u, the 
eigenvalues are obtained from 

h^nUn - fwUn^l = (7) 

With U = {an - a,|^ + /3,|j^) and /.. = (! + |^). Here ql = 2m{V^ - E)/h' 
and ko is the wave vector at z = Zo + 0^. Below, these eigenvalues are used to evaluate the 
corresponding eigenfunction for an specific system. 

I shall now refer to the eigenfunctions \I'^^j,(2;) of finite periodic systems. For this purpose, 
it is convenient first to fix some notation. As shown in figure 1, the coordinates {zj} , with 
j = 0, l,2,...n, define a set of points separated by multiples of Ic, i.e. Zj = jlc- If we 
are interested in evaluating functions at any point z in the j-th cell, it is useful to define 
the difference Sz = z — Zj < Ic- Taking into account the transfer matrix multiplicative 



properties, the total transfer matrix Mt{zo — ^ z) relating the wave vectors ^{zq) and ^{z), 
can be factorized as Mj{z'^ -^ z)Mo^o'{zo -^ z'^), with z'^ = Zo + Sz, and subsequently, the 
wave vector can be written as 



<^(z) 



(l?{z)\ 



( 



aj l3. 



3 t^J 



^ ^ ^(4) ^ 



(8) 



Here aj = Uj — a*Uj^i, [3j = pUj^i, and the functions '0'{z) and Vl-^) represent, depending 
on the difference E — V{z), the right and left propagating or exponentially increasing and 
decreasing functions, respectively [12]. The eigenfunctions are easily obtained by using these 
relations, the boundary conditions and the general expressions for aj and (3j. 

For a scattering system like the one shown in Fig. la), the wave function at z is given 
by 



v^(^,E) = ^«) 






(9) 



It is clear that evaluating the function \E'(z,i?) at E^^^i, we have the corresponding z/-th 
resonant bound state in the /i-th band 



^^^,{z) = ^{z,E, 



1 -^fj^iyj 



(10) 



In figures 3a-d) we plot four different functions for a system with the potential parameters of 
the superlattice GaAs{AlQ^^GaQjAs/GaAsy^. In 3a) we have the function |^2(-2j E)\ for an 
arbitrary energy E (7^ -^^,1^) within the second energy band. This is an extended aperiodic 
wave function with a complicated behavior along the superlattice. In figures 3b and 3c), 
we have instead the functions |^2,2(-2)| and |^E'4,3| corresponding to the second and third 
resonant energies in the second (/i = 2) and fourth (/i = 4) energy bands. In these cases the 
bound-state functions are modulated by an oscillating envelope function with u — 1 minima 
plus those at 2;^ = and Zn = nlc- Notice that, since these wave functions describe not 
only extended but also transmitted states, they are, except for a few points in |^2,2(-2)| , 
different from zero with the probability to find the particle at the two ends of the system 
different from zero. This will change for bounded systems, of course. Finally, in figure 3d) 
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we plot 1^1/ (-2, E)\ at some point in the gap between the second and third bands. The wave 
function behavior, for an incident energy in the gap, is not only compatible with the well 
known vanishing of the transmission coefficient, it makes also evident the localization effect 
induced by the phase coherence at the energy gap, as suggested in some sense by Kohn[2]. 
This is a very appealing result which may deserve further reflections. 

To complete this discussion we shall now briefly refer to the other two types of systems. 
Consider first the system bounded by infinite hard walls. In this case, the /i, z/ eigenf unctions 
are evaluated from 

^^,,{z) = ^(4)(«, +f3*- «* - f3,) (11) 

with E = Efj,^i, and the matrix elements aj and (3j as mentioned above are simple functions 
of the Chebyshev polynomials of order j and j — 1. Care has to be taken when extra pieces of 
length a/2 are added to the system of length nlc- For the particular case shown in figure lb) 
the last two terms in Eq. (11) have to be multiplied by («„ + /9*e"*^'')/(/3„ + Q;*e"*'"'). Two 
of these eigenfunctions (| ^2,2 (-2) I , I ^2,3 (-2) I ) are plotted in figures 4a-b). Contrary to the 
scattering system, the wave function vanishes at the ends. An important difference between 
a periodic potential fianked by infinite hard walls and a constant potential fianked also by 
infinite hard walls, is the presence (in the periodic potential case) or absence (in the single 
well case) of phase interference effect responsible for energy splittings that in the former 
case develop the band structure and the energy eigenvalues. 

Finally, for the system with finite walls mentioned above, the eigenfunctions are given 
by 

^,,,{z) = ^{z',){a,R + P,R + ^{f3,i - a,i)) (12) 

HereEfj^i, is obtained from equation (7). Some of these eigenfunctions are plotted in fig- 
ures 5a-c). Although imperceptible, the wave functions decrease exponentially inside the 
potential walls. As in previous figures, we distinguish two main characteristics: a strongly 
oscillating behavior, which frequency grows with fi and u, and a periodic modulation, with 
period nlc/v. 
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In this letter various new and general formulas for the evaluation of fundamental quanti- 
ties like the resonant energies, bound states, energy eigenvalues and eigenfunctions of finite 
periodic systems, have been reported. Specific and illustrative examples and results, for 
both scattering and bounded one-dimensional systems, have been also presented. 

I acknowledge Project 29026-E from CONACyT, Mexico and the ICTP, Trieste Italy. I 
benefited from comments and discussions with Professors. H. Simanjuntak and N. March. 
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FIGURES 
FIG. 1. Potential profiles, along the direction z, for open, bounded and quasi-bounded n-cell 

systems. 

FIG. 2. Level density in the first band of finite and infinite (Kronig-Penney) 
GaAs (AlQ^^GaQiAs/GaAs)"' superlattice with a = lOOnm, b = 30nm and Vq = 0.23ey. Some 
eigenvalues Ei^u (in eV) for n = 8 are also indicated. 

FIG. 3. Squared wave functions for a system like in Fig. la) but with square barrier potential 
and parameters as in Fig. 2. In a) extended |^^(-z, E)\ with E ^ E^^y and ^u = 2 is shown; In 6), 
and c) the indicated resonant states; In d) a localized wave function in the second gap is shown. 

FIG. 4. Eigenfunctions |^^^jy(z)[ for a system like in Fig. 16) and parameters as in Fig. 2. The 
oscillations frequency and enveloping minima depend on the band and excitation energy indices 
^, V. All functions vanish at the surface. In a) and h) the particle density beneath the surface is 
relatively high. 

FIG. 5. Eigenfunctions in the lower bands for a system like in Figs. Ic) but with square barriers 
and parameters as in Fig. 2), flanked by finite walls. 
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